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Abstract. We compare the results of Balmer-line calculations using recent theory and improved computational 
algorithms with those from the widely-used SYNTHE and BALMER9 routines. The resulting profiles are mostly 
indistinguishable. Good fits to the normalized solar Balmer lines Ha through H<5 are obtained (apart from the 
cores) using the recent unified-broadening calculations by Barklem and his coworkers provided that some ad- 
justment for the continuum is performed. We discuss a surprising linearity with temperature of the Balmer line 
profiles in dwarfs. 

Key words, stars: atmospheres - stars: fundamental parameters 



1. Introduction 

Balmer line strengths are highly sensitive to the temper- 
ature in cool stars because of the 10.2 eV excitation of 
the n = 2 level from which they arise. Fig. 151 from 
Unsold's (1955) classic text illustrates this for H7 equiva- 
lent widths. We show the effect in a different way in Fig. 1, 
based on more recent line-broadening theory. The figure 
is for points on the Ha profile 4 A from the line center, 
but is characteristic of much of the line profile. 

An extensive investigation of Balmer lines in cool 
dwarfs (Fuhrmann, Axer, & Gehren 1993; Fuhrmann, 
Axer, & Gehren 1994) concluded these lines provide 
a more consistent guide to effective temperatures than 
broad-band colors or b — y. Nevertheless, Balmer line pro- 
files are not regularly used to fix the effective temper- 
ature of cool stars. The reasons for this are numerous, 
but have not been explicitly addressed. Some insight may 
be gained from the papers by Van't Veer-Menneret & 
Megessier (1996) or Castelli, Gratton & Kurucz (1997, 
henceforth, CGK). A recent paper which does discuss use 
of Ha in the determination of effective temperatures is by 
Peterson, Dorman & Rood (2001). In addition to the un- 
certainties in placing the continuum level, uncertainties, 
both in the theory of stellar atmospheres (l/H, convec- 
tion) and line formation remain unresolved. 

The absorption coefficient of neutral hydrogen takes 
into account the effects due to the natural absorption 
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Fig. 1. Ha wing strength vs. T e g for several values of log g. 
The profiles are taken from the BP00K2NOVER grid available 
in http://kurucz.harvard.edu 



(natural broadening) , the velocity of the absorbing hydro- 
gen atoms (thermal Doppler and microturbulent broad- 
ening), the interactions with charged perturbers (linear 
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Stark broadening), with neutral perturbers different from 
hydrogen (van der Waals broadening), and with neutral 
hydrogen perturbers (resonance and van der Waals broad- 
ening). Each effect is represented by a profile and the total 
effect requires a convolution. Thermal Doppler and micro- 
turbulent broadenings are described by gaussian functions 
while natural, resonance, and van der Waals broadenings 
have Lorentz profiles. These two profiles are combined into 
a Voigt function. The convolution of the Voigt profile with 
the Stark profile or Stark plus thermal Doppler effect then 
gives the total absorption profile. 

Most of the damping constants and Stark profiles are 
computed from complex theories based on several approx- 
imations, while the complete convolution of all the above 
profiles is a very time consuming algorithm. 

In this paper we describe our attempts to evaluate sev- 
eral aspects of the calculations of Balmer line profiles. 



2. Stark profiles 

Most work on stellar atmospheres makes use of codes pro- 
vided by Kurucz (http://kurucz.harvard.edu). For com- 
puting hydrogen lines the codes are either BALMER9 
(Kurucz, 1993a) which produces profiles for H Q , Hp, H 7 , 
and H<5 or the SYNTHE code (Kurucz, 1993b) which pro- 
duces profiles for any hydrogen line. In the first case Stark 
profiles are interpolated in the Vidal, Cooper, & Smith 
(1973, henceforth VCS) tables, while in the second case 
the Stark profiles are based on the quasi-static Griem the- 
ory with parameters adjusted in such a way that profiles 
from Griem theory fit the VCS profiles of the first mem- 
bers of the Lyman and Balmer series. 

Only the most recent work on the Balmer lines (e.g. 
Barklcm, Piskunov & O'Mara 2000, henceforth, BPO) has 
included the new Stark profiles of Chantal Stehle (hence- 
forth CS) and her coworkers. They are available from a 
link on her website: http://dasgal.obspm.fr/ stehle/. A 
recent reference is Stehle & Hutcheon (1999). 

A problem arises when a given Stark profile is inter- 
polated either in the VCS or in the CS tables by using 
the interpolation method taken from the BALMER9 code. 
This is a bilinear interpolation in log(T) and \og(N e ), fol- 
lowed by a linear interpolation in the parameter Aa = 
A\[A]/F°. Here, F° is the normal field strength in 
Gaussian cgs units, F° = 1.25V?/ 3 , so the interpolation 
in Aa is not independent of the previous one which in- 
volves the electron density N e . We find this introduces a 
small error that shows up as an oscillation in a plot of 
the Stark profile S(Aa) vs depth in the solar atmosphere 
for a small range of displacements from the line center as 
shown in Fig. 2. 

We were able to remove the oscillations by rewriting 
the CS tables with AA as the third (independent) vari- 
able, and using essentially the same interpolation scheme 
as BALMER9. Fortunately, it has resulted that the im- 
proved interpolation leads to no perceptible changes in 
the resulting line profiles. 




Running Depth 

Fig. 2. Normalized Stark width at AA = 0.5A for Ha vs. 137 
depths in an Holweger-Miiller (1974) solar model. Each depth 
step is 0.05 in log(rAsooo)- The vertical lines mark depths cor- 
responding to boundaries of the tables giving S(a) for a fixed 
value of the electron density. 



3. Convolution of profiles and microturbulence 

Neither the BALMER9 code nor the SYNTHE code per- 
form profile convolutions, but all the profiles are simply 
added. In the BALMER9 code, for separations larger than 
0.2 A from the line center, a Lorentz profile (representing 
the natural broadening and the resonance broadening) is 
added linearly to the Stark-thermal Doppler profile inter- 
polated in the VCS tables. For separations smaller than 
0.2 A no Lorentz profile was considered. 

In the SYNTHE code, the Doppler profile, the Stark 
profile, and the Lorentz profile (for natural broadening, 
resonance broadening, and van der Waals broadening from 
He I and H2) are still summed together. The very inner 
core is that of the profile (Doppler, Stark, or Lorentz) with 
the largest full width at half maximum FWHH. 

This method due to Peterson (1993), which we shall 
call the PK approximation, would be rigorously true for 
the wings of two Lorentzians. Since the wing-dependence 
of the Stark profile differs from that of a Lorentzian only 
by ^/(AA), one might expect the approximation to be 
good, as we verified that it is. 

Replacing the sum of the Stark and Lorentz profile in 
BALMER9 by a convolution takes a large amount of com- 
puting time in that the AA step of the convolution has to 
be very small (less than 0.001 A) in order to account for 
the narrow full width at half maximum FWHM of the 
Lorentz profile. This problem can be overcome by includ- 
ing a microturbulent velocity £t in the computations. 

Both the VCS and CS tables include thermal Doppler, 
but not microturbulent broadening. The BALMER9 code 
makes no provision for the inclusion of microturbulence in 
the line profiles owing to the sum of the Stark-thermal- 
Dopplcr profile, interpolated in the VCS tables, with the 
Lorentz profile. The SYNTHE code does allow for a micro- 
turbulence in that it adds the Stark profile to a Doppler- 
microturbulence gaussian profile. 
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Table 1. Models used for Ha tests 



Te(K) 


log g 


^(kms- 1 ) 


Comment 


4500 


1.5 


3.0 


solar abundances 


4760 


1.3 


2.3 


CS22892-052 (cf. Sneden et al. (1996)) 


5770 


4.4 


1.0 


sun 


8000 


3.5 


2.0 


like cool Ap or Am 


8000 


1.5 


12.0 


test of large £ t 


12000 


3.0 


2.0 


hot star 
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Fig. 3. Ha profiles for a model with T eB = 8000^, log 3=1.5. 
The lower curve is for a CSII calculation with an assumed mi- 
croturbulence £t = 12 km s~ . The upper curve, displaced up- 
ward for purposes of illustration, was made using BALMER9, 
the older interpolation scheme for VCS tables, and the PK 
approximation. There is no perceptible difference in the two 
profiles beyond the line core. 



that computed without microturbulence, as is shown in 
Fig. 3. However, the H Q profile computed by SYNTHE 
with no convolutions, but by assuming £t=12 km s _1 
agrees well with the CSII profile. 

The effect of a microturbulent velocity £ t will be small 
until £t approaches the sound speed. It is not surprising, 
therefore, that the only case we have found where plots 
of Ha obtained using BALMER9 with the PK approxi- 
mation and CSII differed significantly is that for £ t of the 
order of the sound speed. Even in this situation, only the 
deepest parts of the core were affected. The line wings still 
matched beautifully. 

The calculations of Fuhrmann ct al. (1993;1994) in- 
cluded Lorentz broadening by a full convolution, while 
BPO used the PK approximation. The above comparisons 
led us to conclude that any differences between their re- 
sults and other calculations (e.g. CGK or Gardiner, Kupka 
& Smalley 1999) cannot be attributed to the PK approx- 
imation or to different Stark profiles (VCS or CS) - the 
immediate line core excepted. 



The only way to rigorously include all broadening 
mechanisms is to do a convolution of the Stark-thermal 
Dopplcr profile, interpolated in the VCS or CS tables, with 
a profile which includes both the Lorentz broadening and 
turbulent motions. If we assume a Gaussian distribution 
of microturbulent velocities, the VCS or CS profiles need 
to be convolved with a Voigt profile. 

To check BALMER9 and SYNTHE profiles we did cal- 
culations using the new CS profiles with improved interpo- 
lation, and a full convolution including a microturbulent 
velocity. We shall refer to such profiles and to the corre- 
sponding code with the abbreviation CSII (Convolution, 
Stehle, improved interpolation). Table 1 shows models pa- 
rameters for which we made calculations of an Ha pro- 
file in order to test the effects of the various approxima- 
tions and improvements mentioned above. All models were 
generated with the ATLAS9 code (Kurucz, 1993a). Solar 
abundances were assumed for all but CS22892-052, for 
which abundances were chosen to roughly match those of 
Sneden et al. (1996). 

We find, with one exception, that the BALMER9 pro- 
files computed with no convolutions and no microturbu- 
lent velocity are in excellent agreement with CSII calcula- 
tions. The only exception occurs for the supersonic micro- 
turbulent velocity £t = 12 km s _1 . In this case the line core 
of the profile computed for £t=12 km s^ 1 is larger than 



4. Broadening of the hydrogen lines by collisions 
with H I atoms 

The BALMER9 and SYNTHE codes allow for the broad- 
ening of the hydrogen lines due to the collisions with other 
neutral H I atoms through the resonance broadening based 
on the Ali & Griem theory (1965,1967). Actually the van 
der Waals effect due to H I should also be included, but 
it can not be simply added to the resonance broadening 
(Lortet & Roueff, 1969) and therefore it was always ne- 
glected in the hydrogen profile calculations. Only recently 
BPO (Barklcm ct al. 2000) presented an unified theory of 
the H I-H I collisions in the stellar atmospheres. The dif- 
ferences in Balmer profiles computed with only resonance 
broadening and with both resonance and van der Waals 
broadenings are fully discussed in BPO. 

We have included in our hydrogen synthetic spec- 
tra (BALMER9, SYNTHE and CSII) the BPO broaden- 
ing. The line half half-width HWHM per unit hydrogen 
atom density w/N{H) is computed according to Anstee 
& O'Mara (1995): 

w/N(H) = (4/7r)"/ 2 r(2 - a/2)va{v Q )(v/v )- a 

where the cross-section a and the velocity parameter a 
for Ha, H/3, and H7 were taken from Table 3 in BPO. 
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Furthermore, we recall that v = (SRT/tt/i) 1 / 2 , where p, is 
the reduced mass for two hydrogen atoms, and w is the 
velocity v for 10 6 cm s _1 . The value of the T function is 
0.901903 for H Q , 0.92437 for and 0.93407 for H 7 . 

In the CSII code, HWHM was computed in accord- 
ing to BPO for each given temperature of the atmo- 
spheric layers. For H,5 the broadening by neutrals was 
obtained by extrapolating BPO's Table 3, but the pro- 
file is dominated by Stark broadening, and is nearly in- 
dependent of the broadening by neutrals. In BALMER9 
and in SYNTHE, HWHM was obtained for each tem- 
perature of the atmospheric layers from a function 
HWHM=HWHM (T/10000) 3 ' where HWHM is the 
value of HWHM for T=10000 K and y was derived from 
the best fit of the above function to the HWHM,T points 
for T ranging from 2000 K to 11500 K at steps of 500 K 
(Fig. 3 in BPO). The parameter y is 0.15 for H„, 0.275 for 
Kp, and 0.30 for H 7 . 

5. Balmer profiles from the Holweger-Miiller solar 
model 

5.1. The solar HM Model 

For the calculation of the solar Balmer profiles we adopted 
the Holweger-Muller model (1974, henceforth, HM) to 
avoid additional complications from various solar mod- 
els, already discussed, for example, by CGK. We started 
from the HM T-r 50 oo relation given for 29 layers, and 
extrapolated-interpolated to suit the depth ranges used 
by our respective codes. 

There are differences in the optical depth coverage of 
the Michigan and Trieste codes. In the first case, the T- 
r 5000 relation was interpolated-extrapolated to 135 layers, 
while in the second case it was interpolated for 50 layers 
before using it in the Kurucz codes. While the Michigan 
code performs integrations directly in terms of log(r5ooo)> 
the use of the Kurucz codes requires a conversion from 
the T5000 depth scale to a RHOX (or / pdx) depth scale, 
where p is the density of the stellar gas and x is the ge- 
ometrical height in the atmosphere. The conversion was 
obtained by computing the continuous opacity K5000 a t 
A=5000 A by means of the ATM code from Holweger, 
Stcffen & Steenbock (1992, private communication) and 
by deriving RHOX from the relation dr 50 oo = K soooP dx. 
The original HM model was made more than a quarter of 
a century ago. Since that time, abundances and the con- 
tinuous opacity routines have been modified, presumably 
for the better. This means that the current relation be- 
tween T5000 and TRosseiand is no longer the same as in the 
HM paper. The latter is inconsistent with the RHOX scale 
of the modern Kurucz codes. 

We adopted as solar abundances the meteoritic values 
from Grevesse & Sauval (1998) and a constant microtur- 
bulent velocity £=1 km s _1 . 

The HM model used in the Kurucz codes is given in 
the Appendix A. 



5.2. Predictions from the HM model 

For clarity, we first list several categories of opacity rele- 
vant to the current problems: 

1. Standard continuous opacity: bound-free and free-free 
transitions in various atoms and ions, Rayleigh and 
Thomson scattering. These are implemented in most 
currently-used model atmosphere and spectrum syn- 
thesis codes. 

2. TOPBASE opacities (Seaton et al. 1992). These opac- 
ities have not yet been widely implemented in cur- 
rent atmosphere codes, so the impact of this important 
work remains to be seen. 

3. Line opacity due to transitions between tabulated 
atomic energy levels. Some of these lines are predicted, 
in the sense that they have not been observed on the 
laboratory, but all relevant levels have been located, 
typically to a fraction of a wavenumber from observed 
lines. We shall call these classified lines. We distinguish 
two categories: 

(a) Stronger lines, which contribute 1% or more to the 
continuous opacity at the central wavelength for 
point in a model atmosphere. 

(b) Weaker lines, for which the above criterion is not 
met. 

4. Line opacity due to transitions involving one and some- 
times two levels whose locations are predicted by an 
atomic structure code. Wavelengths for these lines may 
be uncertain by 10 or more angstroms. A sizable frac- 
tion of these lines involve levels above the first ion- 
ization limit, and the levels are therefore subject to 
autoionization. We shall refer to these as unclassified 
lines. Many of these lines may have been observed in 
laboratory experiments. Again, we list two categories: 

(a) Stronger lines. In certain chemically peculiar stars, 
we know there must be many such lines because we 
are unable to identify a large fraction of the mea- 
surable stellar lines. There are also many uniden- 
tified lines in the solar spectrum, though they are 
usually weaker than a few tens of milliangstroms, 
and typically increase in number to the violet. 

(b) Weaker lines connecting predicted levels. 

5. "Missing" opacity. Calculations of the solar continuum 
using only standard continuous opacity (No. 1 above) 
predict values significantly higher than the "observed" 
continuum. The disparity increases toward the violet 
(see discussion below). 

Fig. 4 compares the solar intensity Ia(0) from the 
center of the Sun measured by Neckel & Labs (1984) 
with Ia(0) predicted using the continuous and line opac- 
ities from Kurucz (1993c) and the HM model given in 
Appendix A. The line opacity is treated with the opacity 
distribution functions (ODF), which include both classi- 
fied and unclassified lines. When such a line opacity is 
considered in the calculations, a rather good agreement 
of the low resolution observations with the low resolution 
predictions is obtained, indicating that much of the miss- 



4 ,. 8 



4 .. 2 



3.8 



* 3.4 



C. R. Cowley and F. Castelli: Balmer lines in the sun and stars 

~i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 





1 ,.4 



1 .. Ol i i i i i i i i i i i i i i i i i i i i i i i i 

3000 3400 3800 4200 4600 5000 5400 5800 6200 6600 7000 7400 7800 

A (A) 



Fig. 4. Comparison of the solar intensity from the center of the sun predicted by the HM model (full line) with the observations 
from Neckel & Labs (1984) (dashed line). The line opacity in this low-resolution calculation is entirely from the ODFs. 



ing opacity could be due to line absorption. Because the 
ODFs involve averages over wavelength intervals of the 
order of 20 A in the 3300-6400 A region and larger for 
A>6400 A, we refer to the calculation of Fig. 4 as a low- 
resolution synthesis. The nature of the missing opacity 
is somewhat controversial, and will not be argued here. 
A recent reference, with citations to earlier discussion, is 
Peterson ct al. (2001). 

Limb darkening predictions from the HM model are 
compared in Fig. 5 with those from Neckel & Labs (1994). 
In this case, opacity from lines is not included in the com- 
putations in accordance with the assumption of Neckel & 
Labs (1994) of observations made at wavelengths free from 
lines contaminating the continuum. The departure of the 
computations from the observations in the violet can be 
explained with the poor chance to have regions free from 
lines in this part of the solar spectrum. Except for the 
violet wavelengths, the agreement is satisfactory. 




Fig. 5. Comparison between observed (points) and com- 
puted (full line) solar limb-darkening curves Ia(cos#)/Ia(0). 
Observations are from Neckel & Labs (1994) and computed 
curves are based on the HM model 



5.3. The Balmer profiles in absolute intensity 

Fig. 6 shows the Balmer profiles for the disk center in abso- 
lute intensity. We have adopted the Kitt Peak observations 
available at the Hamburg site (ftp.hs.uni-hamburg.de; 
pub/outgoing/FTS- Atlas) and described by Neckel (1999, 



henceforth, KPN). The files include absolute intensities, 
as well as continuum estimates at each wavelength. The 
resolution of the observations is about 350000. 

The synthetic Balmer profiles were computed with the 
SYNTHE code and the HM model. Only the relevant 
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Balmer line and standard continuous opacity sources were 
used. The spectrum was degraded at the observed resolu- 
tion and it was broadened by assuming a macroturbulent 
velocity ^ macro — 1.5 km s -1 , although Balmer lines are 
independent of instrumental and macroturbulence broad- 
cnings of the order of those here adopted. 

Fig. 6 indicates that except for H Q , the observed and 
computed profiles can not be directly compared in abso- 
lute intensity, owing to the different levels of the observed 
and computed spectra. The differences are less than 1 % 
for H Q , of the order of 5 % for H^, 4 % for H 7 , and 8 % for 
H^. The discrepancy does not change if all the classified 
lines are included in the computations (see also Fig. 7 in 
CGK). 

The agreement shown by Fig. 4 shortward of 4600 A 
appears better, and apparently contradictory, to that of 
Fig. 6. However, both resolution and opacities are different 
in the two figures. In addition to the averages over the 20 A 
intervals, Fig. 4 includes opacities from unclassified lines. 
When the whole line opacity is averaged over the 20 A 
intervals the detailed spectral differences at each wave- 
length are smoothed off. The result is a better agreement 
of the spectra observed and computed at low resolution 
than that of the same spectra analyzed at high resolution. 

The picture is complicated by absolute measurements 
by Burlov-Vasiljev et al. (1995) of the solar spectral energy 
distribution. It is higher by about 6% than that of Neckel 
& Labs (1984) at H 5 , 4% at H 7 , 2% at R fj , while it is 
about 2% lower at H Q . 

We now turn to a comparison of normalized profiles, 
in which the above problems are less obvious, though nev- 
ertheless present. 

5.4. The normalized Balmer profiles 

In the current work, one of us (CRC) attempted new es- 
timates of the continuum for the observed spectrum - 
less as an attempt to improve on the KPN values, as to 
gain some insight into the uncertainties in this endeavor. 
We began with spectral high points within 10 A intervals 
plotted vs. wavelength, and smoothed the "envelope" by 
selectively deleting points, in an obviously subjective way, 
to achieve an overall smooth plot. The adopted points are 
shown in Table 2, along with those from KPN. We make 
no claim that the current continuum is superior in any 
way to that chosen in KPN. It was simply used in the 
Michigan work for normalization purposes. We employed 
a four-point Lagrange interpolation scheme to normalize 
observations between the chosen points. 

Our independent evaluation of the continuum based 
on the points shown in Table 2 is in excellent agreement 
with KPN, with the exception of the region near H/3. The 
value shown in column 2 for A4861 interpolated with the 
four-point Lagrange formula, from the surrounding points, 
is 1.2% higher than the KPN continuum. This region ap- 
pears depressed for reasons that are unclear and deserve 
investigation. 




M to n« air ' — rtar J — riai 



I [A] 



Fig. 6. Comparison between observed and computed Balmer 
profiles in absolute Ia(0) units. The left scale gives inten- 
sities relative to the following values, all in units of 10 6 
erg cm~ 2 s _1 stear" 1 A -1 . For H Q , the maximum of the up- 
permost panel is 3.2, for H/3 and H7 the maxima are both 5.0, 
and for HS, the maximum of the panel is 5.5. 



The continuous specific intensity using the HM model 
and Michigan codes matches the interpolated continuum 
from Tab. 2 at Ha to within 1%. For H/3 through H<5, the 
calculated continua fall above the measured (as interpo- 
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Table 2. Solar continuum specific intensity in units of 10 cgs 



Wavelength 


This work 


KPN 


3298.973 


0.3235 


0.3231 


3355.431 


0.3269 


0.3272 


3782.919 


0.4083 


0.4093 


4020.705 


0.4589 


0.4591 


4279.262 


0.4652 


0.4666 


4419.404 


0.4598 


0.4609 


4504.079 


0.4540 


0.4545 


4861.000 


0.4230 


0.4179 


5102.095 


0.3999 


0.3990 


5203.252 


0.3906 


0.3902 


5801.460 


0.3435 


0.3424 


6109.561 


0.3200 


0.3189 


6202.178 


0.3146 


0.3144 


6409.847 


0.2990 


0.2972 


6500.584 


0.2907 


0.2899 


6802.324 


0.2660 


0.2663 


6850.076 


0.2619 


0.2627 


6950.356 


0.2546 


0.2553 


6972.875 


0.2536 


0.2540 


7000.000 


0.2524 


0.2524 



Stehle (Stark) + BPO (neutrals) 



lated in Tab. 2) continua by 2.4, 3.9, and 7.8% respec- 
tively. These results agree well with those discussed in the 
previous section of the comparison of the observed and 
computed absolute intensities. 

If we assume the "missing opacity" as cause for these 
disagreements as well as for those shown in Fig. 6, there is 
at present no obviously correct way to account for it. For 
these calculations, we assumed this opacity has the same 
depth dependence as standard continuous opacity sources. 
We have simply scaled them by constant factors until the 
calculated specific continuous intensities agree with the 
observed chosen continuum. 

When spectra normalized to the continuum levels are 
compared, we find an excellent agreement for H Q (Fig. 7). 
The results are the same both from the CSII and the 
SYNTHE code, and are to be compared with BPO's 
Figure 8 (upper), done for the solar flux. We see good 
agreement in all cases. The agreement of the CSII profiles 
with BPO profiles is expected, since the only basic differ- 
ence is the use in BPO of the PK approximation while 
CSII uses a full numerical convolution, a distinction we 
have found thus far to be unimportant. 

As far as the three higher, normalized Balmer lines 
are concerned, the best fits to the wings are obtained 
when the "observed" continua are adjusted downward 
from values obtained by interpolation in Tab. 2 — the sense 
is that the continuum there is too high. For H7 and 
H5, the downward adjustment is 2 %. The observed con- 
tinuum at B./3 needed a downward adjustment of 3 %; 
problems with the continuum in this region were men- 
tioned earlier in this section. Fig. 8 shows the fit for 
H(5. The other two Balmer line fits may be seen at 
the url: http://www.astro.lsa.umich.edu/users/cowley/- 
balmcrs.html/ 
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Fig. 7. Ha profile for the center of the solar disk normalized 
to the continuum level. The thin curve is the observed KPN 
spectrum, and the solid the CSII calculation with an assumed 
microturbulence £t = 1 km s _1 . In this calculation no allowance 
for missing opacity has been made, and the continuum has been 
adopted as described. 
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Fig. 8. KPN spectrum and CSII calculation for H5. 



In principle, the adjustment of the continuum requires 
an iteration with a new continuous opacity to the new 
continuum. Fortunately, the normalized Balmer profiles 
are not very sensitive to small adjustments for the missing 
opacity. 
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6. Inhomogeneities and the Plane Parallel Model 

For perhaps a century we have known that the spectrum of 
the solar photosphere varies from one point on the disk to 
another. The first high-resolution spectra obtained from 
the McMath-Hulbert Observatory showed striking spatial 
variations that came to be known as "wiggley lines." The 
solar line profiles vary markedly, both in time and space, 
and while we have understood the general the nature and 
cause of these variations for decades, recent numerical 
calculations by Nordlund, Stein, and their collaborators 
have provided a detailed description (cf. Nordlund & Stein 
2001). 

In spite of its origin in a turbulent roil, the average line 
spectrum of the sun is remarkably constant. This is partic- 
ularly surprising in the case of the Balmer lines, where the 
large Boltzmann factor (#xiowcr ~ 10) suggests huge local 
non-linear effects. Naively, one would not expect them to 
average out, and the extent to which they do average out 
remains to be fixed. 

In the 1950's, de Jager (1952) attempted to fix the 
temperature fluctuations in the solar atmosphere by mak- 
ing use of the putative nonlinearities of the Balmer lines. 
His conclusions, of temperature differences of a thousand 
degrees from hot to cool columns agrees remarkably with 
modern numerical models. Surely, he was guided by phys- 
ical insight into what the answer needed to be. The Stark- 
broadening theory of that time was rudimentary, and the 
influence of collisions with neutral hydrogen were entirely 
neglected. 

We have found that reasonable matches to the four 
lower Balmer lines can be achieved using modern Stark 
profiles provided recent parameters for broadening by neu- 
tral hydrogen by BPO and the HM model are used. In fact, 
the fits illustrated in Figs. 7 and 8, were all based on the 
empirical plane-parallel Holweger-Muller model, and in- 
clude no attempts to improve the fits by plausible adjust- 
ments of the line-broadening parameters. Other studies 
have explored the sensitivity of the Balmer lines to dif- 
ferent theoretical model atmospheres and to variations in 
the convective mixing length to the pressure scale height 
(l/H). 

We remark here on the surprising linearity of the 
Balmer profiles with the temperature of plane-parallel 
models. This may be illustrated in several ways. In Fig. 1 
we can see that for T e g about 4000K to 6250K the wing 
strengths plot nearly linearly with temperature for the 
three higher gravities. 

This near linearity holds for most points on the line 
profiles, apart from the most central portions. If one takes 
an equally weighted average of Ha fluxes for T = 5500K 
and 6500K, the resulting mean differs imperceptibly from 
that for a 6000K model. Means for a 5000K and 7000K 
model differ only by 2% from the 6000K model beyond 3A 
from the line center. Even for the mean of a 4500K and 
an 7500K model is the difference of the order of 5% (see 
Fig. 9). 
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Fig. 9. Percentage differences in Ha profiles for 6000K model 
and average profiles for three pairs of models as indicated (H a 
profiles from Kurucz 1993a). 



The same effect may be seen in the left panel of Fig. 3 
of Fuhrmann et al. (1993). They show a series of Balmer 
profiles from Ha through US for log <?=4, with effective 
temperatures running from 5000K to 6700K, in steps of 
100K. It can be seen that the different profiles are, for the 
most part, quite evenly spaced. 

The simple means of Fig. 9 are certainly not equiv- 
alent to the detailed calculation performed, for example, 
by Asplund, Nordlund, & Trampedach (1999), based on 
the 3-dimensional numerical models of the solar convec- 
tion zone. Nevertheless, they demonstrate that the non- 
linearities that one might expect from the very large 
Boltzmann factors of the n = 2 level are not realized in the 
resultant Balmer profiles of cool stars. This, in turn, sup- 
ports endeavors to use theoretical profiles from simplified 
stellar models to help fix fundamental stellar parameters. 



7. Conclusions 

We have explored recent techniques for computing Balmer 
line profiles in the sun, and Ha profiles in several mod- 
els with temperatures ranging from 4500K to 12000K. We 
find that new Stark profiles, rigorous convolution, and im- 
proved interpolation techniques make almost no difference 
in the resulting calculated profiles, compared with algo- 
rithms used in the Kurucz codes for several decades. 

Good fits to normalized disk center solar profiles for 
the Ha through H<5 are obtained from the Holweger-Muller 
(HM) model. The Ha profile can also be reasonably fit- 
ted in absolute intensity, but the calculated continua for 
H/3 through H<5 are too high. This may reasonably be at- 
tributed to missing UV opacity, perhaps also to inadequa- 
cies of the HM model used here, as well as to uncertainties 
in the absolute solar calibration. 

In spite of severe temperature inhomogeneities in the 
solar atmosphere, the plane-parallel model appears re- 
markably robust. 
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Appendix A: The Holweger-Miiller (HM) model 
as input for the Kurucz codes 

Table A.f lists the HM model interpolated on 50 depths 
and converted to the RHOX scale of the Kurucz codes. 
The last four columns are the input model for the Kurucz 

codes. The units for the continuous opacity K5000 arc m 
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Table A.l. The HM solar model 
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